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Abstract 

We investigate near threshold neutral pion photoproduction off protons to fourth order in heavy 
baryon chiral perturbation theory in the light of the new data from MAMI. We show that the 
unitarity cusp at the secondary 7r+n threshold is in agreement with expectations from the hnal 
state theorem. We also analyze the fourth order corrections to the P-wave low-energy theorems 
and show that potentially large A-isobar contributions are cancelled by sizeable pion loop effects. 
This solidifies the parameter free third order predictions, which are in good agreement with the 
data. 



1 Introduction 



Chiral perturbation theory is the tool to systematically investigate the consequences of the sponta- 
neous and explicit chiral symmetry breaking in QCD. S-matrix elements and transition currents of 
quark operators are calculated with the help of an effective field theory formulated in terms of the 
asymptotically observed fields, the Goldstone bosons and the low-lying baryons. A systematic per- 
turbative expansion in terms of small external momenta and meson masses is possible. We call this 
double expansion from here on chiral expansion and denote the small parameters collectively by q. Be- 
yond leading order, coupling constants not fixed by chiral symmetry appear, the so-called low-energy 
constants (LECs). These have to be determined by a fit to some data or using some model. A large 
variety of processes such as pion-nucleon scattering, real and virtual Compton scattering and so on 
has already been investigated in this framework, sharpening our understanding of the chiral structure 
of QCD (for reviews, see e.g. ||]). 

Neutral pion photoproduction off protons and deuterons (which gives access to the elementary neutron 
amplitudes) is one of the prime processes to test our understanding of the chiral pion-nucleon dynamics 
for essentially two reasons. First, over the last decade fairly precise differential and total cross section 
data have been obtained at MAMI |^, ^, ^ and SAL |6|, ^] . A further experiment involving linearly 
polarized photons was performed at MAMI, which not only improved the differential cross sections 
but also gave the first determination of the photon asymmetry Q. Second, the S-wave amplitude i?o+ 
is sensitive to a particular pion loop effect [l^. In the threshold region, the fourth order heavy baryon 
chiral perturbation theory calculation (HBCHPT) (which involves the sum ai + a2 of two low-energy 
constants) agrees with what is found in the multipole analysis of the data js], 11, 12|. In addition. 



the rather counterintuitive prediction for the electric dipole amplitude for photoproduction off the 



^0 



neutron, |-E'o+"i > 1-^0+ I translates into a threshold deuteron amplitude |13[ that was verified by 



^0+ I ^ 1-^0+ 

a SAL experiment within 20% Moreover, in [11| it was also shown that there are two P-wave low 
energy theorems (LETs) for the Pi^2 multipoles which show a rapid convergence based on the third 
order calculation. While the LET for Pi could be tested and verified from the unpolarized data, only 
the recent MAMI measurement of 7p ^ vr^p allows to disentangle the contribution from the P2 and 
the P3 multipoles (the latter being largely determined by the LEG hp at third order). ^ It has been 
frequently argued that contributions from the delta isobar, that only appear at fourth order in the 
chiral expansion for Pi and will not only spoil the rapid convergence of the P-wave LETs but 
also lead to numerically different values. This is witnessed e.g. in an effective field theory approach 
including the delta as an active degree of freedom p5[| in which one counts the nucleon-delta mass 
splitting as another small parameter. A third order analysis in that framework seems to indicate 
indeed large corrections rendering the agreement of the prediction for Pi at threshold with the value 



deduced from the differential cross sections as accidental 116 



In this paper, we complete the fourth order (complete one loop) analysis based on HBCHPT by 
evaluating the corresponding corrections for the three P-wave multipoles. We use this framework 
to analyze the new data from MAMI, which confirms and sharpens previous findings concerning the 
electric dipole amplitude Eqj^ and sheds new light on the convergence issue of the P-wave LETs. 



*^Note that the first but somewhat model-dependent comparison between P-wave multipoles and the LET predictions 
was given by Bergstrom p^ . 
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2 Formal aspects 



In this section, we collect some basic formulas needed for describing the reaction j{k) + p{pi) — > 
7^^{q) + p{p2)- In the threshold region, it is legitimate to consider -k^ photoproduction in S- and 
P-wave approximation, with the corresponding transition current matrix element given by 

T ■ e = ia ■ e [Eo+{lo) + k ■ q Pi{uj)] + ia ■ k e- q P2{uj) + {q X k) ■ e Psiu) . (1) 



Here, m = 938.27 MeV is the proton mass, s = {pi + k)"^ = {p2 + q)^ the total centre-of-mass 
(cm) energy squared , a; = (s — m? + M^o)/2^/s the cm energy of the produced neutral pion, and 
€^ = (0, e ) the polarization vector of the real photon in the Coulomb gauge subject to the transversality 
condition e • A: = 0. At threshold, the vr^ is produced at rest in the cm frame, g = 0, so that 
(Jthr = M^o = 134.97 MeV corresponding to \/stiir ~ + The secondary 7r"'"n-threshold opens 
at (jJc = 140.11 MeV where y/s^ = M^++mn (with m„ = 939.57 MeV the neutron mass). At that point, 
the strong unitary cusp related to the rescattering process 7p ^ vr'^'n ^ ir^p occurs in the electric 
dipole amplitude Eqj^{u). In the vicinity of the cusp its generic form reads Eq-\-{uj) = —a— 6^/1 — uJ^ joJ^ 
with two constants a and h. The amplitudes -Pi, 2, 3(1^) are linear combinations of the more commonly 
used magnetic dipole {M\±) and electric quadrupole (i^'i+) P-wave pion photoproduction multipoles. 
The combinations -Pi, 2, 3(1^) arise most naturally from the decomposition of the T-matrix in eq.(l). 
Of importance for the later discussion are also the threshold P-wave slopes Pi, 2, defined via 

Pi,2 = lim — -- — , (2) 

<?-*o \q\ 



because for these the LETs mentioned in the introduction have been derived in |llj]. The differential 
cross section and the photon asymmetry S can be expressed in terms of the multipoles as 



I A: I da 



\q \ dflc 



A + B cos9 + C cos^ 9 , (3) 



A = \Eo+\'+'^\P2\' + lm' , (4) 
B = 2Re(So+A*) , (5) 
C = |Pi|2-i|P2p-^|P3|' , (6) 



2\k\ \d^lcm 

with 9 the cm scattering angle and we have dropped the argument lo. 



3 Chiral expansion of the multipoles 

We wish to calculate the T-matrix element, eq.(l), to order 0{q'^). For the electric dipole amplitude 
£'o+(w), this has been already done in |11]. In that paper, also the third order terms for the P-wave 
multipoles Pi,2,3(w) were evaluated. Here, we give these up-to-and-including fourth order. We make 
use of the standard heavy baryon effective chiral Lagrangian, which has been given to complete one 
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loop accuracy, i.e. 0{q^), in ref.[^. The expressions for the multipoles spht in three parts. First, one 
has the (renormahzed) Born terms which subsume the lowest order couplings {g-KN^i^) complemented 
by the anomalous magnetic moment (Kp) contributions together with all pion-loop corrections of 
these parameters at order 0{q^) and 0{q^). Secondly, there are the pion loop graphs with at most 

(2) 

one insertion from the dimension two Lagrangian C)^^- For these, the one-nucleon reducible parts 
which just renormalize the Born terms are taken out. Thirdly, there are the one-nucleon irreducible 
counterterm contributions which lead to simple polynomial amplitudes. 

Consider first the renormahzed Born terms which are expressed in terms of the physical parameters 
Qt^]^ ,m, Hp. The second and third order terms for Pi2^^{uj) are given in the appendix of ref.|ll]. We 
display here only the novel fourth order contributions. 



jBorn 



(^) 



3207rm4 \ 



(21 + 19Kp)M; + (20 + QKp)oj' 



(43 + lbKp)Ml - (26 + 10Kp)a;2| , 



>Born I 



3207rm4 I a;2 

^{(3.,-2)M.^-(13 + 8.>^} 



(8) 

(9) 
(10) 



leOvrm^ 

with \q \ = ^Juj'^ — M'^. From now on = 134.97 MeV denotes the neutral pion mass, g.,^N = 13.1 
is the strong pion-nucleon coupling constant and Kp = 1.793 the anomalous magnetic moment of the 
proton. 

Next, we give the P-wave contributions from the fourth order one loop graphs. According to their 
prefactor, gA or g\, these fall into two classes. In these loop diagrams charged as well as neutral 
pions occur in internal lines and we have neglected throughout the small mass difference — 
M^o = 4.6 MeV. In sharp contrast to the S-wave amplitude £'o+ (having a strong cusp effect) this 
approximation is legitimate for P-wave amplitudes since their imaginary parts and consequently their 
cusp effects are extremely tiny corrections. The numerical differences which result from taking the 
charged or neutral pion mass for M^^ should be regarded as an intrinsic inaccuracy of our 0{q'^) 
calculation. 

First, we give the analytical expressions for the C'((?^)-loops proportional to gA'- 
p!°°P(c^) 



egA \q I 
m(47rF^)3 

8C4 



{ (3 + ^C4)a;2 In ^ + \u^MI - + ^ arcsin^ ^ 



3uj 



w^)^/^ arcsin 



UJ 



jloop 



(c) 



ddA \q\ 



+- 



TT 



ijj 



Ml -u? - Ml arcsin 



M^ 
A 

UJ 



UJ 



2uj-^ 



oj^ arcsm — - H == 

^M^ 



UJ 



UJ 

- arccos 

? M^ 



(11) 



^(4vrF.)3 1 + r^y 1- T + + + ^^'^ - 2^^^^ ™' 



M^ 



4ca 



2uj + ^{uj' + 2M'^] 
6u) 



Ml 



uj^ arcsin ■ 



UJ 



jloop 



(^) 



egA \q\ 



IT 



(1 + 4C4) 



UJ 



Ml-uj^ + Ml arcsin 



UJ 
'M~n 



(12) 
(13) 
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with = 92.4 MeV the weak pion decay constant, Kn = —1.913 the neutron magnetic moment, 
QA = g-KNF-K/mn = 1.29 and C4 = mc/i. The low-energy constant C4 has been determined from pion- 
nucleon scattering inside the Mandelstam triangle as C4 = 3.4 GeV^^ |l^. A is the scale of dimensional 
regularization which will be set equal to A = m. As a check these loop contributions fulfill the condition 
p|°2 3(0) = which confirms that all one-nucleon reducible pieces are indeed taken out. The analytical 
continuation above threshold u > is obtained by the following substitutions: 

UJ 



Ml 



Ml. 



arcsm — > h « In — — 

M^ 2 M^ 



Ml 



(14) 



Similar analytical expressions are found for the other class of 0((7^)-loops proportional to 



jloop 



egj \q\ 
vrM2 



— (7 + 4Kp + 2«;n)ln— ^ 



3u 



(7 + 4Kp + 2Kn) arcsin — — 



+ ■ 



6^3 
2 



4M^ - Guj^M^ 



^ + 3u:{2u:' - Ml) arcsin ^ + {lu^ - AM^^M^-oj^^ 



+^(17 + 8Kp + 4k„) - y(19 + WKp + 5Kn) + ^(^' - Ml) arcsin^ ^| , (15) 



eg A \q\ 
m(47rF^ 



^(5 + 2Kp + Kn)uj'^ In ^ + ^^(19 + WKp + 5k„ 
o A 9 



+ 



2,uj 



Ml{7 + iKp + 2k„) - 2uj^{5 + 2Kp + k„) 



arcsin 



^_Ml 
M^ 2cj2 



UJ 



arcsm 



+ ■ 



ttM^ 
e^i 1^1 



4M^ - + 3w^ - 4M^(m2 



a; 



2\3/2 
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11 + 8Kp + AKr, 



m(47rF^)3 {3lo 



IT 



— (^Kji 3Kp 3) 



Mt 



{Ml - u;2)3/2 



M^ 
(16) 

(17) 



which also fulfill the nontrivial condition Pi°2^{0) = 0. 

Finally, we are left with the polynomial counterterm contributions, 



Pf(^) 
P2\u;) 



egA \q 1 1^ 
m(47rF^)3 

egA \q\ ^'^ 
m{4TTFT,)^ 



6(A), 



6(A), (18) 

(19) 

e|,|6p{.-^}. (20) 

The introduced new parameters (LECs) 6,2 (A) are dimensionless and they balance of course the 
scale dependence appearing in the fourth order loop contribution via the chiral logarithm ln(M7r/A). 
The form of P^^{ijj) follows from the relativistic operators Og and O9 constructed in ref.|jl^. The 
LEG hp already appeared at third order, the only new feature here is a kinematical correction |A; | = 
oj — Ml 127X1 + . . ., which at threshold amounts to a 7% reduction. 

Furthermore, we give the resonance contributions to the low energy constants ^1^2 and hp. As mentioned 
in ref.|l^] there is a small contribution to bp from t-channel vector meson exchange (/9'^(770) and 
01(782)) and consequently also an analogous vector meson exchange contribution to ^1^2, 

(47rF.)3' ^1 - - 



(21) 
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Here, we have used various simplifying relations (see ref.|Tl|) for the vector meson coupling constants 
together with the KSFR relation for the vector meson masses Mp ~ M^. The dominant contribution 
to hp and ^1^2 come from the low-lying A(1232) resonance. In refs.dn], ||] we used a relativistic 
tree level approach in which the delta contribution is parametrized in terms of four couplings gi, 
§2, Y and Z. The latter two are so-called off-shell parameters emerging in a relativistic description 
of the spin-3/2 fields. In a corresponding effective Lagrangian, these would be represented by some 
higher order contact interactions. In fact, most of the delta resonance physics can be represented 
by the static isobar approach, which can be thought of as the leading term in a systematic effective 

it was 



A 



field theory expansion like the one given in |15|. In the study of pion-nucleon scattering |H 
already demonstrated that the dominant isobar contributions come indeed from the lowest order Born 
graphs. Therefore, we use here the static non-relativistic isobar-model where one gets the following 
expressions, 

, (A) _ K* QA 



(A) 



"?2 



(22) 



(23) 



3^2 A2 - M2 

with A = 293 MeV the delta-nucleon mass-splitting and n* the A^A transition magnetic moment. It is 
important for the numerical evaluation to keep the M^-term in the denominator, this is also justified 
in the small scale expansion, where one counts A as a small parameter like the pion mass. 

For the later discussion of the P-wave LETs, we now give the various contributions to the slopes -Pi,2- 
We start with the renormalized Born terms expressed by the physical pion-nucleon coupling constant 
g-KN, the renormalized anomalous magnetic moment Hp and the proton mass m. Furthermore, we 
introduce the small parameter ^ = M^o/m ~ 0.144 and have 

,•2 

1+ Kr- 



Pi(Born) 
P2(Born) 



egjrN 



Svrm^ 



|(2 + Kp) + ^(9 + 5«p) 



I - Hp + -{?, + Kp) - —(13 + bKp) 
Z o 



(24) 
(25) 



Here, novel terms of order /x^ appear. The result for the chiral loops at order 0{q^) can be taken from 



Pi((?3-loop) 

P2(g^-l00p) 



3847r2m2 



(10 - 37r) 



1927r2m2 



(26) 
(27) 



These contributions are known to be quite small. From the formulae for p|°2^(l(j) and Pf 2(1 
above, one can readily deduce the terms due to the chiral loops and counterterms at order 0{q' 



LO) given 

4\ 



Pi (g"^— loop, ct) 



m{i7rFT,) 



8 o2 

-C4 + 3 + ^(7 + 4Kp + 2K„) 



In- 



2 ^ 2 

+(1 + 2g\)^ + C4 + 2 + (13 + 4Kp + 2k, 
o y io 



A 

) + ei(A) 



2 5vr 

gAY 



P2(g'' — loop, ct) 



m(47rP^) 



:C4 



vr 



-(4c4 + 2 + 5l)— + 



2 ■ 
40 



25i 



(5 -|- t^n) 



In 



vr 



C4 + 1 + ^(5 
9 18^ 



4k^ 



2K„) + e2(A) 



(28) 



(29) 
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The resulting numerical values will be given later. 



4 Results and discussion 



The new MAMI differential cross section data span the energy range from = 145.1 MeV to 
165.6 MeV in steps of about 1.1 MeV. In addition, the photon asymmetry S has been evaluated 
for energies from 145 to 166 MeV, all these data for S have been binned to one average energy of 
= 159.5 MeV. In total, we have 171 differential cross section data, 19 total cross section points and 
7 data points for the photon asymmetry S. 

First, is it instructive to compare the new data with the previously obtained ones of Fuchs et al. Q. 
For doing that, we compare fits using the fourth order expressions for the S-wave amplitude Eq^ and 



the third order ones for the P-wave amplitudes (as it was done in |11, 12|). The resulting LECs and 
X^/doi are collected in table ffl. In both cases the two S~wave LECs are completely anticorrelated. 





Schmidt et al. 


Fuchs et al. 


ai [GeV~^] 


10.585 


3.464 


a2 [GeV-4] 


-4.542 


3.136 


Corr(ai, 02) 


-0.998 


-0.999 


ai + 02 [GeV-4] 


6.04 


6.60 


bp [GeV-3] 


14.84 


13.00 


No. of data 


190 


180 


xVdof 


3.19 


2.20 



Table 1: Values of the LECs resulting from a three parameter fit to 
the cross section data of ref. Q and [^. Corr. denotes correlation 
between the two S-wave LECs. 



i.e. only the sum ai + 02 is of relevance. It agrees within 10% for the two fits, showing that the 
S-wave multipole £'0+ is internally consistent. The P-wave LEG bp is somewhat increased, but now 
consistent with the value obtained from fitting the SAL data |^, ~ 15GeV~^ [|^. This can 

simply be traced back to the fact that the new MAMI total cross sections are larger than the old ones 
above the secondary yr^n-threshold. It is gratifying that this so far puzzling experimental discrepancy 
is now resolved. 

Next, we wish to investigate the strength of the S-wave cusp. For that, we use the realistic two- 
parameter model developed in ref.[|ll| (which is similar to the so-called unitary fit of ref.|^]), where 
i?o+ is given by 

Eo+ {uj) = -a-b y^l- ^ . (30) 

Assuming isospin invariance for vrN-rescattering, the strength of the cusp given by the parameter 
h = \pla~ -E'o+" can be inferred from the well measured pion-nucleon scattering length a{7r~p 
ir^n) and the precise CHPT prediction for the electric dipole amplitude i?o+(7J' ~^ 7r~^n) at threshold 
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Figure 1: The electric dipole amplitude in the threshold region. Left panel: Real part. The 
various lines are explained in the text. Right panel: The modulus of i?o+ for the full fit (see 
text) in comparison to the SAL data |^ 0. 



(which agrees with the data). This gives b = (3.67 ± 0.14) • 10~^/iV4+ [§. Fitting the older MAMI 
data, the resulting value for b came out sizeably smaller, b ~ 2.8- 10~^ /M^+ [12|. This prompted some 



speculations that the strength of the unitary cusp is very sensitive to isospin violation. If, however, 
we use this same model together with the third order predictions for the P-waves and apply it to the 
new MAMI data, we get 



a = 0.54 • IQ-^/M^^ 
b = 3.63 • 10"V^7^^ 



b„ 



14.43 GeV" 



(31) 



with a x^/dof of 3.21, which is of the same quality as the one of the three parameter HBCHPT fit 
discussed before. The value for b in eq.(^) is in perfect agreement with the prediction obtained from 
the final state theorem and assuming isospin invariance for vrN-rescattering. That sheds some doubt 
on the speculation that a precise measurement of the unitary cusp would be a good tool to investigate 
isospin violation. The resulting real part of £'o+ is given by the dash-dotted line in fig. ||. In this 
simple three parameter approach, one can also predict the photon asymmetry S, since P3 is governed 
by the LEG bp and P2 is given by the LET. Of course, in the full fit involving also the fourth order 
contribution to the P-waves, one has additional terms from loops and counterterms, nevertheless, in 
this simplified ansatz one can already estimate the corrections to be expected from these additional 
terms. Using the LECs collected in table |l], we obtain the dash-dotted in line fig. |2|, which agrees quite 
well with the data from the MAMI analysis. Therefore, we conclude that the corrections to the P2 
multipole should be small. Also, we had already noted before that P3 is only modified on the percent 
level by the fourth order corrections. Thus, to keep the fine balance between and |i-*3p, which 
governs the size of S, only modest corrections to P2 should be expected. 
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Figure 2: The photon asymmetry at = 159.5 MeV. The 
various hues are explained in the text. The data are from |^]. 





Set I 


Set II 


ai [GeV-^] 


7.734 


8.588 


as [GeV-4] 


-1.506 


-2.288 


Corr(ai, 02) 


-0.998 


-0.998 


ai + 02 [GeV-4] 


6.23 


6.30 


xVdof 


1.36 


1.35 


^0+(^thr) [W~yM^+] 


-1.13 


-1.12 


Eo+{lOc) [10'VA4+] 


-0.53 


-0.52 



Table 2: Values of the S~wave LECs and £'0+ at the two thresh- 
olds resulting from the five parameter fits of the data of ref.Q. 
Corr (01,02) denotes correlation between the two S-wave LECs. 

We now discuss the full fits including the fourth order corrections to the P-wave multipoles Pi^2,3(^)- 
We have performed two types of fits. In set I, we only fit to the differential and total cross section 
data excluding the photon asymmetry S. For set II, we include the photon asymmetry data in the 
fits. Let us first discuss the electric dipole amplitude. It should come out (largely) independent 
of the fitting procedure since S is only indirectly sensitive to the S-wave. The resulting LECs are 
collected in table As expected, one finds a very similar result for ai + 02 in agreement with the 
previous determinations. The resulting £'o+('^) comes out indeed independent of the fitting procedure 
as shown by the solid and dashed lines in fig. 0. It is in good agreement with determinations based 
on the older MAMI and the SAL data, which lead to £;o+(t^thr) = -(1-3 ± 0.2) • 10"^/A4+ and also 
with the result obtained from the new MAMI data, Eo+iaJthr) = -(1.33 ± 0.08 ± 0.03) • 10"^/M^+. 
Therefore, even though the convergence of the chiral expansion in this multipole is slow, the fourth 
order calculation is able to describe it in the threshold region with one parameter (the sum of LECs 
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01+02)- This small value for £'0+ at threshold clearly establishes the large pion loop effect first pointed 
out in [|lO|. For completeness, we also show in fig. |l| the modulus of the electric dipole amplitude, 
I £'0+1 = ([Re £0+]^ + [Im -£'0+]^)^''^) in comparison to the data from SAL |^], which nicely shows 
the unitary cusp. 

We now turn to the P-waves. Here, we encounter the following problem. While the best fit of type I 
gives a good x^/dof, see table there is an almost perfect correlation between bp, ^1 and ^2 (which 
is expected since the differential cross sections are only sensitive to + a-nd the resulting 

values for bp or ^2 come out either too large (based on expectations from resonance exchange, to be 
discussed below) or with too large uncertainty. We have also performed fits with fixing 6p at the 
previously determined value of 14.8 GeV~^, which gives almost the same x^/dof but a wastly different 
value for ^2- Furthermore, including the leading effects of D-waves in the low energy region does not 
change this. On the other hand, in all cases the LEG ^1 comes out in a narrow range, which is in 
agreement with the estimates based on resonance exchange to be discussed next. If one includes the 
photon asymmetry data, the cross sections and the photon asymmetry are well described, see the solid 
line in fig. ^, but the resulting value for bp is too large whereas ^1 and .^2 come out of the size expected 
from resonance saturation. For completeness, we show the energy dependence of the real parts of the 
P-wave multipoles for the fit including the asymmetry data in fig. ^. 

8.0 
6.0 
4.0 

^ 2.0 



DC 

-4.0 
-6.0 
-8.0 

0.140 0.145 0.150 0.155 0.160 0.165 

EJGeV] 

Figure 3: Real part of the P~waves Pi (i = 1, 2, 3) versus photon energy. 




Therefore, to get a more reliable estimate of the corrections to the P-wave LETs, we employ the 
resonance saturation hypothesis. First, taking the parameters used here, the LET predictions read 
(which are nothing but the sum of the third order renormalized Born and loop terms): 

P^^^ = 0.469 GeV~2 , [0.445, 0.492] GeV^^ , 

P^^^ = -0.498 GeV-2 , [-0.472, -0.523] GeV'^ , 
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where the numbers in the square brackets refer to a 5% theoretical uncertainty 
the new MAMI data are [g] 

P^^P = (0.441 ± 0.004 ± 0.013) GeV'^ , 
P^^P = (-0.440 ± 0.005 ± 0.013) GeV-^ , 

which are in good agreement with the LET predictions. From our fourth order results, we get for the 
sum of renormalized Born, third and fourth order loop and counterterm contributions 

Pi = (0.460 + 0.017 - 0.133 + 0.0048^1) GeV"2 , (34) 
P2 = -(0.449 + 0.058 - 0.109 + 0.0048^2) GeV"2 , (35) 

where the ^1^2 only depend on the A^A transition magnetic moment. We note the rather sizeable (25%) 
correction from the fourth order loops which at first sight seems to destroy the agreement between the 
LETs and the data. However, it is known that k* ~ 4 ... 6, so we collect in table ^ the predictions for 
Pi^2 for reasonable variations of k* . We see that for k* = 4, the delta contribution almost completely 
cancels the large fourth order loop effect and thus the predictions for the P-wave slopes are within 
7% of the empirical values. Note, however, that the empirical finding |Pi^^| = 1^2^^ | is difficult to 
reconcile with any theory. 



K* 


Pi [GeV-2] 


P2 [GeV-2] 


4.0 


0.408 


-0.475 


4.5 


0.416 


-0.487 


5.0 


0.427 


-0.498 


5.5 


0.439 


-0.509 


6.0 


0.450 


-0.521 


6.5 


0.461 


-0.532 



Table 3: Prediction for the P-wave slopes for varia- 
tions of the A A transition magnetic moment n* . 



. The results based on 
(33) 



5 Summary 

In this paper, we have studied near threshold neutral pion photoproduction off protons in the frame- 
work of heavy baryon chiral perturbation theory to complete one loop (fourth order) accuracy, updating 
and extending previous work on this topic |11, 12|. The pertinent results of this investigation can be 
summarized as follows: 

(i) We have given the fourth order corrections (loops and counterterms) to the three P-wave mul- 
tipoles Pi, 2, 3- Two new low-energy constants appear, one for Pi and the other for P2. We have 
also given analytic expressions for the corrections to the low-energy theorems for the P-wave 
slopes Pi, 2 5 see eqs.(28,29). 
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(ii) We have analyzed the new MAMI data first in the same approximation as it was done in 
previous works (i.e. the P-waves to third order only). Using a realistic two parameter model for 
the energy dependence of the electric dipole amplitude Eq+, we have extracted the strength of 
the unitary cusp which agrees with the prediction based on the final state theorem. 

(iii) Using the full one loop amplitude, one can fit the cross section data and the photon asymmetry. 
The combination of S-wave LECs is stable and agrees with previous determinations, leading to 
-E'o+('^thr) = —1-1 • 10^'^/M^+. Two of the three P-wave LECs are not well determined because 
of strong correlations. More photon asymmetry data are needed to cure this situation. 

(iv) We have analyzed the new LECs in the framework of resonance saturation in terms of (dominant) 
A-isobar and (small) vector meson excitations. The isobar contributions depend only on the 
strength of the A^A transition magnetic moment. We have shown that for reasonable values of 
this constant, the 25% fourth order loop corrections to the P-wave LETs are almost completely 
cancelled by the isobar terms. This solidifies the third order LET predictions, which are in good 
agreement with the data, cf. eqs. 
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